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ABSTRACT 

We report the discovery of TrES-4, a hot Jupiter that transits the star GSC 02620-00648 every 
3.55 days. From high-resolution spectroscopy of the star we estimate a stellar effective temperature of 
Toff- — 6100 ± 150 K, and from high-precision z and B photometry of the transit we constrain the ratio 
of the semi-major axis a and the stellar radius i?* to be a/i?* — 6.03 ±0.13. We compare these values 
to model stellar isochrones to constrain the stellar mass to be A/^ = 1.22 ± 0.17 M©. Based on this 
estimate and the photometric time series, we constrain the stellar radius to be i?^, = 1.738 ± 0.092 Rq 
and the planet radius to be Rp = 1.674 ± 0.094 i?jup. We model our radial-velocity data assuming a 
circular orbit and find a planetary mass of 0.84 ± 0.10 Mjup. Our radial- velocity observations rule out 
line-bisector variations that would indicate a specious detection resulting from a blend of an eclipsing 
binary system. TrES-4 has the largest radius and lowest density of any of the known transiting planets. 
It presents a challenge to current models of the physical structure of hot Jupiters, and indicates that 
the diversity of physical properties amongst the members of this class of exoplanets has yet to be fully 
explored. 

Subject headings: planetary systems — techniques: photometric — techniques: radial velocities — 
techniques: spectroscopic 



1. INTRODUCTION 

Despite the ever increasing number of discovered tran- 
siting exoplanets (18 at the time of writing), our un- 
derstanding of the relationships between host star prop- 
erties and the planets' physical and orbital parame- 
ters is still incomplete. While the mass-radius rela- 
tion for most transiting pla nets agrees with the models 
(jCharbonneau et al.ll2007"al ). several planets have either 
too high or too low d ensities. In particular, HD209458b 
(iKnutson et all [200l. HAT-P-lb (f Bakos ct alj |2007j) , 
and WASP-lb (|Collier Cameron et al. 2007) all have 
much larger radii and lower densities than expected 
for planets of their mass and distance from the host 
star. Several mechanisms for explaining this dis- 
crepancy have been propo sed, among them differ- 
ent internal heat sources (|Guillot fc ShowmanI [200l 
iBodenhei mer et all l2003f) . incr eased planetary atmo- 
spheric opacities ()Burro-ws et al.li2007 ). and less efficient 
heat tr ansport as a result of int erior composition gra- 
dients (jChabrier fc Baraffd l2007l ). Our Trans-atlantic 
Exoplanet Survey (TrES) has previously announced the 
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discovery of three transiting planets (jAlonso et al.ll200l : 
lO 'Donovan et all2006l[2b07 ). each with distinctive prop- 
erties. We describe here the newly discovered transit- 
ing planet TrES-4, whose mean density of p = 0.222 ± 
0.045 g cm"-^ is the lowest of all known exoplanets with 
measured radii and masses. 

2. PHOTOMETRY AND SPECTROSCOPY 

We monitored a 5?8 x 5?8 field in Hercules with 
the Low ell Observatory Plane t Search Survey Telescope 
(PSST, iDunham et "all l2004l ) and the Sleuth telescope 
at Palomar Observatory between UT 2006 May 6 and 

2006 August 2. All images were processed and the pho- 
tometry and t r ansit search carried out as described in 
iDunham et al.l (|2004D . Both telescopes detected transits 
of the host star GSC 02620-00648: PSST observed 2 full 
and 1 partial transits, and Sleuth observed 3 full and 
4 partial transits. The shape of the events was consis- 
tent with the transit of a Jupiter-sized planet across an 
F dwarf, and we undertook a program of foUowup obser- 
vations to confirm the planetary nature of the object. 

We observe d the candida te with the CfA Digital 
Speedometer (jLathaml |1992[ ) from 2006 September to 

2007 April. We obtained 7 spectra covering 45 A cen- 
tered at 5187 A, with a resolving power of A/AA ^ 35 000 
and S/N between 11 and 13 per resolution element. This 
spectral region includes the gravity-sensitive Mg b triplet 
and the host star properties derived from an analysis 
of this region will depend on the star's metallicity. We 
obtained the radial velocities (RVs) by cross-correlation 
against a synthetic template chosen from a large li- 
brary of calculated spectra based on Kurucz model atmo- 
spheres (sec Nordstrom ct al. 1994; Latham et al. 2002). 
These velocities have a typical precision of 0.5 kms~^ 
and they show no significant variation within the er- 
rors. We derived the effective temperature (Tes) and 
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projected rotational velocity (wsini) by comparing our 
observed spectra against synthetic spectra w ith a wide 
range of parameters (see [Torres et al.l |2002|) . The val- 
ues of Tcff and v sin i are listed in Table [1] and assume 
[Fe/H] = 0.0. We also estimate the surface gravity 
to be log g — 3.8 ± 0.2. We combined the Johnson- 
Cousins photometry that we obtained (see below) with 
archive 2MASS data, a nd used the color-temperatur e cal- 
ibrations for dwar f s by|Ramfrez fc Melended (l2005h and 
ICasagrande et al.| {2006) to estimate Teg. The average 
result, Tes = 6130 ±80 K, is consistent with the spectro- 
scopic value. 

The mass and radius of the star were determined on 
the basis of the spectroscopic T^s, the value of a/i?* de- 
rived from the light curve fit described below, and an 
assumed metallicity of [Fe/H] — 0.0 ± 0.2. The quan- 
tity a/ J?* is closely related to the stellar density, and is 
determined in this case with much higher relative pre- 
cision than l ogff. It is therefore a better proxy for lu- 
minosity (see Sozz etti et ahl 120071 ). We searched for the 
best match between stellar evolution models interpolated 
to a fine grid in age and metallicity and the three ob- 
servables, within the stated errors. Uncertainties were 
derived from the full range in A/^ and Ri, allowed by 
the isochrones within observational err ors. A compari - 
son with the models from the series bv lYi et all (l20?)lh 
yielded = 1.22 ± 0.17 Mg, = 1.738 ± 0.092 Rq, 
and logg = 4.045 ± 0.034. This value of \ogg is consis- 
tent with, but better constrained than the spectroscopic 
estimate. The age we derive is 4.7 ± 2.0 Gyr, and the 
predicted absolute visual magnitude {My — 3.36 ±0.27) 
implies a distance of 440±60 pc, ignoring extinction. The 
error estimates above do not include possible systematic 
errors in the stellar evolution models. The gravity, radius 
and age estimates for the host star indicate that this is 
an evolved star at the base of the subgiant branch. 

In order to characterize the host star, we obtained 
off-transit BV{RI)c photometry of TrES-4 on UT 2007 
April 14 with the 1.05-m Hall telescope at Lowell Obser- 
vatory in combination with a 2K x 2K SITe CCD. We 
calibrated the photometry using 7 standard fields from 
iLandoltl (|1992D . The results are hsted in Tablc[T]together 
with other relevant data for the host star of TrES-4. 

We carried out high-precision in-transit z-band pho- 
tometry of TrES-4 oi l UT 2007 M ay 3 and UT 2007 May 
10 with KeplerCam (|Holman et al. 2006) at the Fred L. 
Whipple Observatory (FLWO) 1.2-m telescope, and B- 
band photometry on UT 2007 May 10 using NASACam 
at the Lowell Observatory 0.8-m telescope. With Ke- 
plerCam we gathered 373 and 540 z-band 30-second ex- 
posures on UT 2007 May 3 and UT 2007 May 10, re- 
spectively. With NASACam, we obtained 192 60-second 
S-band exposures. For both data sets, we derived differ- 
ential fluxes relative to an ensemble of local comparison 
stars. The photometry is shown in Fig. [TJ 

We obtained high-precision RV measurements of TrES- 
4 on UT 2007 March 27-29, usi ng HIRES and its I2 
absorption cell (jVogt et al.l Il994l ) on the Keck I tele- 
scope. Our spectra have a nominal resolving power 
A/AA ~ 55 000 and a typical signal-to-noise ratio S/N ~ 
120 pixel^^. We used nine echelle orders in the wave- 
length range 3200-8800 A to derive the velocities for 
TrES-4. The four star-plus-l2 spectra (plus one l2-free 
template exposure) provide good coverage of the critical 



TABLE 1 
TrES-4 Host Star 



Parameter 


Units 


Value 


Source 


RA 


J2000.0 


17''53™13!05 


1 


Decl. 


J2000.0 


+37°12'42'.'6 


1 


GSC 




02620-00648 






mas yr^-'^ 


[-6.5,-23.5] 


1 


V 




11.592 ± 0.004 


2 


B - V 




0.520 ± 0.007 


2 


V - Rc 




0.312 ± 0.008 


2 


V-Ic 




0.601 ±0.008 


2 


J 




10.583 ± 0.018 


3 


J-H 




0.233 ±0.024 


3 


J-Ks 




0.253 ±0.026 


3 




Mq 


1.22 ±0.17 


2 


R*^ 


Rq 


1.738 ±0.092 


2 




Lq 


3.74 ±0.86 


2 




K 


6100 ± 150 


2 


logg 


cm s~^ 


4.045 ± 0.034 


2 


vsini 


kms"'^ 


9.5 ± 1.0 


2 


Age 


Gyr 


4.7 ±2.0 


2 


Distance 


pc 


440 ± 60 


2 


References. — (1) 


UCAC2 (Zacharias ct al. 



[200i ): (2) this paper; (3) 2MASS iSkrutskie et al. 



^ This value of ij* is derived from fitting the en- 
tire light curve and is better constrained than the 
value obtained solely from the stellar evolution mod- 
els (Sec. 2). The uncertainty in i?^ includes both the 
statistical incertainty of 0.044 Rq and the 14% un- 
certainty in Mir 
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Fig. 1. — High-precision foUowup 2-band and B-band photome- 
try of TrES-4. The plot shows the relative flux (including the color- 
dependent extinction correction) of the TrES-4 system as a function 
of time relative to the center of transit, adopting the ephemeris in 
Table [3] Each light curve is labeled with the telescope and date of 
observation. The residuals from the simultaneous fits (solid lines) 
are shown below each light curve. 

phases. As the ephemeris of the system is fixed, RV data 
are needed only to determine the amplitud e and systemic 
veloc ity of the orbit. In these cases (see iKonacki et ahl 
|2003() a few measurements are enough. We extracted and 
reduced all raw spectra using the A4AKEE software writ- 
ten by T. Barlow fsee iSozzetti et al. 2006 1. The final RV 
values are listed in Table [2l The typical precision of the 
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TABLE 2 
Radial Velocity 
Measurements of TrES-4 



40 r 

20 - 

r 
-20 r 
-40 r 



0.0 0-2 0.4 0.6 0.8 1.0 

Phase 

Fig. 2. — Top: Radial velocity observations of TrES-4 obtained 
with Keck/HIRES using the I2 cell, shown relative to the center 
of mass and adopting the ephemeris in Table |3] The best-fit orbit 
(solid line) is overplotted. Middle: Residuals from the best-fit 
model to the radial velocities. Bottom: Bisector spans shifted to a 
median of zero, for each of the iodine exposures as well as for the 
template (which is shown as an additional point at phase 0.669). 

radial velocities {(Jrv) — 10 ms^^ is slightly degraded 
by the modest stellar rotation. The high-resolution, high 
S/N Keck spectra were then used to rule out blend sce- 
narios (see below), and they are presently being analyzed 
for improved characterization of the host star. 

We fit a Keplerian orbit to these data assuming zero 
eccentricity as a good first approximation, as expected 
from theoretical arguments for a period as short as 
3.55 days. The period and epoch of transit were held 
fixed. The rms of this fit is 11 ms^^, which is similar 
to the internal errors of the velocities. The parameters 
of this orbital solution are listed in Table [3] We find 
Mp sini = 0.731 ± 0.057 [{M^ + Mp)/Mq]2/3 Mjup. The 
orbit is displayed in Fig. [2] (top panel) along with the ob- 
servations, and residuals are shown in the middle panel. 

We investigated the possibility that the RV variations 
are the result of distortions in the line profiles caused 
by contamination f rom an unresolved eclipsing binary 
(jSantos et al.ll2002l: [Torres et alll2005D . instead of being 
caused by a planetary companion. We cross-correlated 
each Keck spectrum against a synthetic template match- 
ing the properties of the star, and averaged the corre- 
lation functions over all orders blueward of the region 
affected by the I2 lines. From this representation of the 
average spectral line profile we computed the mean bi- 
sectors, and as a measure of the line asymmetry we cal- 
culated the "bisector spans" as the velocity difference 
between points selected near the t op and bottom of the 
mean bisectors ([Torres et al.ll2005l ). If the RV variations 
were the result of a blend with an eclipsing binary, we 
would expect the line bisectors to vary in phase with the 
photometric perio d with an amplitude s imilar t o that 
of th e velocities (jQueloz et al.l [20011; (M andushev et al.l 
|2005[ ). Instead, we detect no variation in excess of the 
measurement uncertainties (see Fig. [21 bottom panel), 
and we conclude that the RV variations are real and that 
the star is orbited by a Jovian planet. 

3. PROPERTIES OF TRES-4 AND DISCUSSION 



HJD 


RV (m 


s-l) 


2454187.07793 


97.5 ± 


11.8 


2454188.11122 


58.9 ± 


9.4 


2454189.02631 


-76.5 ± 


11.6 


2454189.14920 


-73.9 ± 


10.3 



We analyzed the two z-band and one S-band pho- 
tometric time series u sing the analytic light curves of 
iMandel fc Agol ([2002| ). We assumed a circular orbit and 
a quadratic stellar limb-darkening law, fixing the coeffi- 
cients at the color-dependent values tabulated in iClaretl 
(|2000l '20041 for the spectroscopic T^s and logg, and as- 
suming solar metallicity. We first estimated the time of 
center of transit Tc by fitting a model light curve to the z 
data gathered on UT 2007 May 10. We then determined 
the orbital period P by combining these data with the 
TrES discovery data (which affords a baseline of 1.0 yrs). 
We then fixed the values of Tc and P (listed in Table [3]) 
in the subsequent analysis. 

Our model has six free parameters: the planet radius 
i?p, the stellar radius i?*, the orbital inclination i, and 
the color-dependent extinction for each of the three pho- 
tometric datasets, kzi, kz2, and ks- We assume that the 
observed flux is proportional to e^'"'", where m denotes 
the air mass. The values of Rp and R^^ as constrained 
by the light curves alone are covariant with M^. In our 
analysis, we first estimated the quantity a/R^, which is 
independent of the assumed value of M^, and then used 
this estimate to constrain the va l ue of M^, from ste l- 
lar isochrones ([Sozzetti et al.ll2007l : [Holman et al.ll2007D . 
We then fixed = 1.22 Mq, and estimated the sys- 
tematic error in the radii using the scaling relations 
1/3 

i?p cx i?^ oc Mj ( see footnote to Table [T]). 

We found first the values of Rp, i?^, z, k^i, kz2, and 
ks that minimized the using the AMOEBA algo- 
rithm (Press ct al. 1992). This model is shown as the 
solid curves in Fig. [1] We then conducted a Markov 
Chain Monte Carlo analysis similar to tha t describ ed in 
Holman et all (I2006D . iCharbonneau et all ([2007b[ ). and 



Winn et all (|2007[ ). We created two MCMC chains with 



400,000 points each, one starting from the best-fit values 
and one starting from a random perturbation to those 
values. We then rejected the first 23% of the points to 
minimize the impact of the initial conditions and found 
the results from the two chains to be indistinguishable. 
We examined the histograms of the six parameters, as 
well as the histograms for several combinations of pa- 
rameters relevant to anticipated foUowup studies. We 
assigned the optimal value to the median, and the 1-a 
errors to the symmetric range about the median than en- 
compassed 68.3% of the values. We list these estimates 
in Table [1 

TrES-4 has the largest radius and the lowest density 
of all exoplanets whose mass and radius are known, and 
as such presents n ew challenges f or the theory of irra- 
diated gas giants. iBurrows et al.l (120071 ) suggested that 
increased planetary atmospheric opacities and the inclu- 
sion of a transit rad ius correction (Baraffe et al.l l2003t 
IBurrows et al.l l2003f ) can explain the radii of all large 



4 



TABLE 3 
TrES-4 Planet Parameters 





Units 


Value 


P 


days 


3.553945 ± 0.000075 




HJD 


2 454 230.9053 ± 0.0005 


a 


AU 


0.0488 ± 0.0022 


i 


deg 


82.81 ±0.33 


a/Ri, 




6.026 ±0.131 


b = a cos i 1 Ri, 




0.755 ±0.015 


K 


ms~^ 


97.4 ±7.2 


'Shires 


m s"'^ 


+23.7 ±5.8 


Mp 


Mjup 


0.84 ±0.10 


-ftp 


R}\ip 


1.674 ±0.094 


P 


g cm~^ 


0.222 ± 0.045 


logs 


cm s"'^ 


2.871 ±0.038 


Rp/ R* 




0.09903 ± 0.00088 



" The uncertainty in Rp includes both the statistical 
incertainty of 0.053 iJ,jup and the 14% uncertainty in 
Mi, from Table [1] 

planets. It appears, however, that TrES-4's radius is 
stiU too large for its mass, age and insolation. In terms 
of mass and distance to the host star TrES-4 is similar to 
HD209458b, but has nearly 30% larger radius (i?p = 1.67 
vs. R-p — 1.32). Because of the higher host star luminos- 
ity TrES-4 receives about twice the stellar flux, but it is 
unlikely that the radius difference can be explained solely 
by the higher stellar irradiation. For example, HAT-P- 
Ib, which has about the same radius as IID209458b and 
is the only other planet with p < 0.3 gcm~^, receives 
only 60% of the flux that HD209458 b does. None of the 
recent models of ho t Jupiters {e.g., iFortnev et al.|[2007t 
iBurrows et "alll2007[ ) can predict a radius as large as that 
of TrES-4 at its orbital separation for the estimated age 
and for any mass, even when higher atmospheric opaci- 
ties are considered and allowance for the transit radius 
correction is made. 

The properties of TrES-4 and its host star allow for 
many interesting foUowup studies, some of which are al- 
ready underway. The large radius of TrES-4 in the visual 



suggests that it may have an extended outer atmosphere, 
sir nilar to that detected arou nd IID209458b in Lyman a 
bv lVidal-Madiar et al.l ()2003D . If such an extended enve- 
lope is caused by mass loss through the planet's Roche 
lobe boundary, TrES-4 may have an even bigger enve- 
lope because of its smaller Roche limit: i?R, nche ^ 3.5 Rr, 



vs. 



R 



Roche 



4.4 Rp for HD209458b (jErkaev et all 



|2007() . Moreover, TrES-4 might be in a strong hydro- 
dynamic "blow-off' regime where the outer atmospheric 
layers are detached from the planet's gr avitationa.1 field 
and escape in a comet-like tail (iVidal-Madiar et al.ll2003t 
iLecavelier des Etangs et al.ll2Q04[ ). 

The relatively fast rotation {vsmi = 9.5kms~^) and 
brightness of the TrES-4's host star, as well as the 
planet's size, are favorable for measuring the Rossiter- 
McLaughlin effect (RME) {e.g., ClaudiJjJWinn 2007). 
Of particular interest is the angle between the star's spin 
axis and the planet's orbital plane, as it can provide in- 
formation about the possible mechanisms of the planet 
migration and interaction with the protoplanetary disk. 
The semi-amplitude of the RME is proportional to v sin i 
and to {Rp/Ri:Y and we estimate that for TrES-4 the 
RME can reach about 90 ms""'^, or roughly equal to the 
semi-amplitude of its orbital velocity. 
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